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Objective: Estimate the potential impact of Avahan, the India AIDS Initiative, among 
female sex workers (FSWs) and their clients in five districts of Karnataka state, south India. 

Design: Examination of time trends in sexually transmitted infection (STI)/HIV preva- 
lence from serial cross-sectional surveys, combined with mathematical modelling. 

Methods: Survey data from each district were used to monitor changes in FSW STI/HIV 
prevalence during Avahan. A deterministic model, parameterized with district-specific 
survey data, was used to simulate HIV/HSV-2/syphilis transmission among high-risk 
groups in each district. Latin hypercube sampling was used to obtain multiple parameter 
sets that reproduced district-specific HIV prevalence trends. A Bayesian framework 
tested whether self-reported increases in consistent condom use (CCU) during Avahan 
were more compatible with FSW HIV prevalence trends than assuming no or slow 
(preintervention rates) CCU increases, and were used to estimate HIV incidence and 
infections averted. 

Results: Declines in FSW HIV prevalence occurred over 5 years in all districts, and 
were statistically significant in three. Self-reported increases in CCU were more 
consistent with observed declines in HIV prevalence in three districts. In all five 
districts, an estimated 25-64% (32-70%) HIV infections were averted among FSWs 
(clients) over 5 years. This corresponded to 142-2092 FSW infections averted depend- 
ing on the district (two-fold to nine-fold more among clients). 

Conclusion: Empirical HIV prevalence trends combined with Bayesian modelling have 
provided plausible evidence that Avahan has reduced HIV transmission among FSWs 
and their clients. If current CCU levels are sustained, FSW HIV prevalence could decline 
to low levels by 2015, with many more infections averted. 
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Introduction 



The HIV epidemic in southern India is highly 
concentrated among high-risk groups (HRGs), particu- 
larly female sex workers (FSWs), their clients and MSM, 
who contribute disproportionately to HIV transmission 
compared to the general population [1—3]. Given this 
epidemiological context, Avahan, the India AIDS 
Initiative, was developed by the Bill & Melinda Gates 
Foundation and implemented from December 2003 
[2-4]. Avahan's objective was to scale up targeted 
interventions to at least 80% of HRGs in Karnataka, 
Andhra Pradesh, Maharashtra and Tamil Nadu states, to 
reduce HIV transmission among them and subsequently 
in the general population [1—4]. The common preven- 
tion packages included peer-led outreach and education; 
free condoms; sexually transmitted infection (STI) 
treatment services; referrals for HIV and tuberculosis 
testing, and HIV care; empowerment and community 
mobilization to address local structural barriers; and mass 
communication programmes. Other services, including 
condom and STI treatment social marketing and 
behaviour change communication, were also targeted 
specifically to clients of FSWs [3—7]. 

As previously explained [3,8,9], the main Avahan impact 
evaluation on behavioural and HIV outcomes relies on 
Integrated Behavioural and Biological Assessments 
(IBBAs surveys) among targeted HRGs in selected 
intervention districts, complemented with mathematical 
modelling. 

This paper evaluates the potential impact of Avahan on 
HIV transmission among FSWs and clients in five districts 
of Karnataka State - Bangalore urban, Bellary, Belgaum, 
Mysore and Shimoga, with different epidemiology and 
history of intervention. Karnataka state has an estimated 
76 000 urban FSWs, of whom 60 000 are covered by 
Avahan in 18 of 27 districts (Supplementary material, 
Box li, http://links.lww.com/QAD/A319) [10,11]. In 
Mysore and Shimoga, Avahan was and remains the first 
and only intervention implemented among HRGs 
[3,10,11]. In Belgaum and Bellary, Avahan was not the 
first HRG intervention but is now the only one. It is the 
major HRG intervention but not the only one in 
Bangalore Urban. Shimoga, Bellary and Bangalore Urban 
districts have lower FSW HIV prevalences compared to 
Belgaum and Mysore [6,7]. In these five districts, 
intervention activities were initiated and scaled up rapidly 
between January 2004 and June 2005, which prevented 
true baseline measurements [7,9,12,13]. As no compar- 
able HIV data among FSWs in Karnataka are available for 
non-Avahan districts, there is no comparable control 
group [3,7,8,13-15]. 

The following main analyses are undertaken for each 
district. First, we determine whether decreases in FSW 
STI/HIV prevalence have been observed following 



Avahan implementation. Second, HIV/STI transmission 
modelling is calibrated (using behavioural and prevalence 
data) and used within a Bayesian framework to determine 
whether observed HIV prevalence trends support the 
evidence of self-reported increases in condom use during 
Avahan (hypothesis testing) and, if so, to estimate the 
likely impact of those increases in condom use on HIV 
prevalence, incidence and infections averted among 
FSWs and clients. The modelling results and strength 
of evidence provided by the hypothesis testing are 
discussed and interpreted in the light of additional 
evidence on implementation, coverage and intensity of 
the intervention and relationship between exposure of 
FSWs to the programme and associated changes in 
condom use. 



Methods 

Time trend analysis 

HIV/STI trend analysis is based on serial rounds of 
IBBAs, involving random samples of FSWs, using 
traditional and time location cluster sampling [6,7]. 
Three rounds of IBBA have been conducted in all five 
districts between 2004 and the end of 201 1 (Table 1). The 
first IBBA (Rl) was conducted 7—19 months after 
programme initiation. Round 2 (R2) and Round 3 (R3) 
follow-up surveys occurred 28—37 and 56—73 months 
after Rl, respectively. Crude and adjusted (for variables 
that differed significantly between rounds in each district) 
logistic regression analyses were performed using STATA, 
version 10.0 (Stata Corp., College Station, Texas, USA) 
and took into account the sampling weights and cluster 
sampling design. Logistic regression was used, instead of 
binomial regression, to avoid convergence issues when 
controlling for many confounding variables. The study 
was approved by the institutional review boards of St 
John's Medical College in Bangalore, India, and the 
University of Manitoba in Winnipeg, Canada. Additional 
details on data collection and statistical analyses are 
provided in Supplementary material (Methods A), http:// 
links.lww.com/QAD/A319 [6,7], 

HIV/sexually transmitted infection mathematical 
model 

Model description 

A deterministic compartmental model was used to 
simulate HIV/herpes simplex virus type 2 (HSV-2)/ 
syphilis (Tp) transmission between FSWs and their clients 
through commercial sex and through longer-term 
noncommercial partnerships [13]. As chlamydia (R2: 
2.8—12.0%) and gonorrhoea (R2: 1.3—3.1%) prevalence 
was deemed low, urine samples were not collected at R3 
and, thus, these infections were not modelled (Table 1). 
FSWs were stratified by duration of selling sex (overall 
duration varied between 4 and 14 years across districts), 
two levels of commercial sexual activity (low and high) 
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and three levels of consistent condom use (CCU) with 
occasional clients ('never', 'sometimes' and 'every time'), 
which varied over time (Supplementary Table SI, http:/ / 
links.lww.com/QAD/A319). The fraction of sex acts 
protected by condoms was specified for each CCU level; 
low and high activity FSWs had specific frequencies of 
commercial sex, which could be varied across IBBA 
rounds. An adjustment factor, estimated using comp- 
lementary data, was also applied to account for over- 
reporting of FSW condom use (Table 2, Supplementary 
material, Methods B, http://links.lww.com/QAD/ 
A3 19) [16J. Condom efficacy was also taken into account 
[13,17—22]. Occasional clients were stratified by the 
duration and frequency of commercial sex. Given the lack 
of sexual mixing data between FSWs and clients, random 
mixing was assumed according to availability (i.e. 
proportionately to total number of commercial sex 
partnerships reported by each risk group). 

HIV was modelled with an initial phase of high- 
infectivity, a long low-infectivity phase and a pre-AIDS 
phase of increased infectivity [23—26]. Those with AIDS 
were assumed to cease being sexually active. Cocirculat- 
ing HSV-2 was modelled dynamically with associated 
cofactors representing facilitation of HIV and HSV-2 
acquisition and transmission [27—29]. Access to anti- 
retroviral treatment (ART) was limited until recently (~8 
and 20% of eligible patients with CD4 < 250 by end of 
2005 and 2007, respectively); hence, it was not modelled 
[30,31]. Syphilis was modelled dynamically, including a 
transient immune stage and disease progression, with 
background treatment before Avahan, and Avahan- 
specific treatment and periodic presumptive treatment 
(PPT) [32-34]. Other than the stage of HIV and HSV-2/ 
Tp cofactors, the force of HIV/STI infection also 
depended on gender, frequency of sex acts with different 
partnership types and condom use (details of the model 
and parameters in Table 2, supplementary material 
(Methods B), http://links.lww.com/QAD/A319, and 
supplementary tables S1-S2, http://links.lww.com/ 
QAD/A319). The mathematical model was used to 
simulate STI/HIV transmission in each district as follows. 

Parameter assumptions 

Plausible uniform ranges for the prior distributions of 
biological and district-specific behavioural parameters 
were specified based on literature reviews and the FSW 
and client IBBA data, respectively (Table 2, Supple- 
mentary tables S1-S2, http://links.lww.com/QAD/ 
A319). 

Given the lack of baseline and historical measurements, 
the model's prior distribution for time trends in CCU (i.e. 
fraction of FSWs who are 'every time' users) during 
commercial sex with occasional clients between the start 
of the epidemic and the start of Avahan was based on 
estimated CCU trends by Lowndes et al. [9] - 
reconstructed by utilizing two IBBA survey questions 



that asked FSWs when they started using condoms 
consistently and when they started commercial sex 
work [9,13] (Supplementary material, Boxlii, http:// 
links.lww.com/QAD/A319). Following the start of 
Avahan, CCU with occasional clients in presence of 
the intervention (i.e. 'estimated CCU trends') was 
assumed to increase linearly up to levels reported in 
each IBBA survey and to remain constant after the last 
round. The validity of these trends is addressed in the 
discussion. The proportion of FSWs 'sometimes' using 
condoms after the start of the intervention was assumed to 
increase linearly from 0% to the estimates from each IBBA 
round (<10% in each district). The proportion of FSWs 
'never' using condoms varied to ensure the sum of 
proportions always equals 1 . In absence of the interven- 
tion, as no empirical control groups were available, two 
plausible control groups (i.e. counterfactuals) were 
defined for trends in CCU covering the period following 
Avahan initiation. 'Control 1' assumed no change in 
CCU (fixed at the CCU level at the start of Avahan), 
whereas 'Control 2' (more conservative than Control 1) 
assumed CCU would have continued to increase slowly 
at the preintervention rate suggested by the 'estimated 
CCU trends' analysis (Supplementary figure S5, http:// 
links.lww.com/QAD/A319). Based on available data, 
condom use in noncommercial long-term relationships 
between HRGs was assumed to be low and constant over 
time. Syphilis treatment and PPT were parameterized 
using Avahan 's programme monitoring data. 

Fitting process and plan of analysis 
Latin hypercube sampling was used to randomly select 
multiple parameter sets from the prior parameter 
distributions. The model was run with each parameter 
set between 1987 and 2015. Posterior parameter sets were 
selected if a set produced modelled FSW and client HIV 
prevalence projections that lay within the 95% confidence 
interval (95% CI) of the FSW and client HIV prevalence 
in all IBBA rounds and prespecffied FSW HSV-2/Tp 
prevalence ranges (for R2 and R3, HIV prevalence 
figures were adjusted on the Rl distribution of 
prespecified key confounding variables, excluding those 
specifically varied over time in the mathematical model) 
(Supplementary material, Methods C, http://links.lww. 
com/QAD/A319). 

To test whether self-reported increases in condom use 
during Avahan were consistent with observed HIV 
prevalence trends (hypothesis testing), the fitting process 
was repeated for each CCU trend hypothesis (Estimated, 
Control 1, Control 2). Following previously used 
methods [13], the 'estimated CCU trend' hypothesis 
was deemed most likely only if it achieved a higher 
frequency of fits (proportion of prior parameter sets that 
produced model fits) than the two alternative control 
CCU trend hypotheses. The hypothesis testing stages 
were undertaken to validate the CCU trends as they were 
retrospectively estimated from self-reported CCU data. 
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Table 2. Ranges for district-specific model parameters (prior distribution) for condom use and syphilis treatment before and after Avahan 
implementation sampled from during the model fitting stage. d 



Definition of model input 



Mysore 



Belgaum 



Bellary 



Shimoga 



Bangalore 



Condom use between FSWs and clients for each condom use time trends hypothesis before and after the start of Avahan (Supplementary figure S5, http://links.lww.com/ 
QAD/A319) 



'Estimated CCU trends' (in presence of intervention) 



Annual rate of increase in condom use 


prior to intervention 


0.016- 


-0.038 


0.033- 


-0.057 


0.039- 


-0.065 


0.027-0.060 


0-0.044 


Proportion who are consistent condom 


users 


















Start of the intervention 




0.13 


-0.24 


0.42- 


-0.55 


0.30- 


-0.48 


0.21-0.39 


0.39-0.54 


IBBA R1 




0.31 


-0.40 


0.87- 


-0.95 


0.67- 


-0.79 


0.49-0.62 


0.74-0.82 


IBBA R2 




0.62 


-0.73 


0.86- 


-0.93 


0.85- 


-0.93 


0.79-0.89 


0.83-0.89 


IBBA R3 




0.88- 


-0.94 


0.95- 


-0.99 


0.83- 


-0.92 


na 


na 


Control group 1 




















Annual rate of increase in condom use 


prior to intervention 


0.016- 


-0.038 


0.033- 


-0.057 


0.039- 


-0.065 


0.027-0.060 


0-0.436 


Proportion who are consistent condom 


users from the start of the 


0.13- 


-0.24 


0.42- 


-0.55 


0.30- 


-0.48 


0.21-0.39 


0.39-0.54 



intervention onwards 
Control group 2 

Annual rate of increase in condom use prior to intervention until 

the latest IBBA (R3 except in Shimoga and Bangalore) 
Proportion who are consistent condom users at the start of 

the intervention 
For all condom assumptions (CCU, Control 1 and 2) 

Proportion of FSW who report 'sometimes/often' using condoms 

with occasional clients 

IBBA R2 

IBBA R3 

Proportion of sex acts with occasional clients when a condom is used by 

FSW 'Always' using condoms 

FSW 'Often/sometimes' using condoms 

FSW 'Never' using condoms 
Adjustment factor for overreporting of condom use - 'relative reduction 

in the fraction of commercial sex acts protected by condoms b 
Fraction (%) of sex acts protected in main partnerships between 

married or cohabiting partners before and after 2004 
Condom efficacy per sex act for HIV (%) [1 7,1 8] 
Condom efficacy per sex act for HSV-2 and syphilis (%) [13,29,30] 
Syphilis treatment 

Fraction (%) of FSW receiving antibiotics in absence of Avahan 

(i.e. background treatment) 
Fraction {%) of clients symptomatic for syphilis c 
Fraction (%) of clients seeking background treatment when 

symptomatic 

Fraction {%) of FSWs seeking Avahan treatment for syphilis 3 
Number of FSWs receiving presumptive positive treatment (PPT) a,c 

From start of Avahan to Jan 2007 (number receiving 

PPT every 3 months) 

After Jan 2007 (number receiving PPT every 6 months) 
Mean duration of infection of FSW seeking Avahan treatment 

for syphilis (months) 
Assumptions for impact assessment 
Start of Avahan 

Consistent condom use in presence of intervention 

Condom use trends (counterfactual) in simulated control group 

after 2004 in absence of intervention 
Syphilis treatment in control group 



0.016-0.038 0.033-0.057 0.039-0.065 0.027-0.060 



0.13-0.24 



0.42-0.55 



0.30-0.48 



0.21-0.39 



January 2004 April 2004 July 2004 July 2004 

Based on district-specific 'estimated CCU trends' 
Control 1 or 2 



0-0.436 



0.39-0.54 



0.25- 


-0.37 


0.07- 


-0.14 


0.06- 


-0.13 


0.07-0.17 


0.09-0.15 


0.06- 


-0.12 


0.01- 


-0.05 


0.05- 


-0.11 


na 


na 


0.82- 


-0.93 


0.92- 


-0.97 


0.91- 


-0.97 


0.83-0.94 


0.91-0.97 


0.54- 


-0.67 


0.34- 


-0.84 


0.44- 


-0.77 


0.53-0.74 


0.56-0.74 


0.07- 


-0.38 


0.0- 


-0.63 


0.34- 


-0.77 


0.18-0.46 


0.02-0.32 


0.00- 


-0.25 


0.00- 


-0.25 


0.00- 


-0.25 


0.00-0.25 


0.00-0.25 


2.0- 


-13.9 


4.7- 


-13.3 


3.0- 


-9.4 


2.0-6.4 


8.6-16.6 


80- 


-95 


80- 


-95 


80- 


-95 


80-95 


80-95 


40- 


-70 


40- 


-70 


40- 


-70 


40-70 


40-70 


9.8- 


-27.8 


9.8- 


-27.8 


9.8- 


-27.8 


9.8-27.8 


9.8-27.8 


21- 


-50 


21- 


-50 


21- 


-50 


21-50 


21-50 


0- 


-100 


0- 


-100 


0- 


-100 


0-100 


0-100 


10.0- 


-30.0 


10.0- 


-30.0 


10.0- 


-30.0 


10.0-30.0 


10.0-30.0 


132- 


-175 


282- 


-431 


979- 


-1501 


80-163 


281-810 


435- 


-633 


595- 


-646 


117- 


-1021 


13-89 


1000-1656 


0.26- 


-0.40 


0.46- 


-0.71 


0.58- 


-0.81 


0.40-0.65 


0.36-0.54 



June 2005 



No Avahan PPT after 2004, 0% seeking Avahan Tp treatment 



Additional model parameters and fitting data shown in Supplementary tables S1 -S3, http://links.lww.com/QAD/A31 9. 'na', not available at the time of the modelling 
analysis - R3 data were only used in trends analysis (Table 1), which are quicker to carry. CCU, consistent condom use; FSW, female sex worker; IBBA, Integrated 
Behavioural and Biological Assessment. 

a Avahan programme monitoring (MIS) data; and calculation as in [21]. 

b The adjustment factor is based on the ratio of condom use reported by clients vs. female sex workers in IBBA [16]. 
c Lower limit based on IBBA data across districts, upper limit based on Orroth et a/. [22]. 



Ranges are generally based on 95% confidence intervals derived from the IBBA surveys. 
e Reynold Washington personal communication. 



Our hypothesis testing method was validated on 
simulated data and showed good discriminating power 
(Supplementary material: Methods D, http://links.lww- 
com/QAD/A319). The posterior parameter sets from 
the 'estimated CCU trends' fitting results were then used 
to simulate the trends in HIV prevalence and incidence in 
presence of intervention, and in its absence by substitut- 
ing the Control 1 or 2 CCU trends (all other parameters 



remaining the same). The number and prevented fraction 
of HIV infections since the start of Avahan were estimated 
from these simulated incidence trends with multiple 
posterior parameter sets being used to reflect the 
uncertainty in the model projections, using 95% 
credibility intervals (95% Crls). The impact estimates 
using Control 2 are more conservative than those using 
Control 1 because Control 2 assumes CCU increases in 
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the absence of an intervention. With either control group, 
no additional syphilis treatment by Avahan was assumed. 



Results 

Data analysis of time trends 

HIV prevalence declines occurred across rounds in all 
districts but were only statistically significant over 5 years 
(Rl to R3) in Mysore, Belgaum and Bellary in the 
adjusted analysis (Table 1). Chlamydia, gonorrhoea and 
syphilis prevalence also declined significantly between Rl 
and R2 or R3 in Mysore. Adjusted STI trends, between 
Rl and R2, were more variable in other districts. 
Declines in chlamydia prevalence were significant in 
Shimoga and Belgaum, nonsignificant in Bellary, but 
increased in Bangalore. Gonorrhoea declined signifi- 
cantly in Belgaum and Bellary but not significantly in 
Shimoga or Bangalore. No significant trends in syphilis 
prevalence occurred except in Mysore. 

Mathematical modelling 

Hypothesis testing of trends in condom use 
The 'estimated CCU trends' was most likely for Bellary, 
Shimoga, Mysore and Bangalore, with only Bangalore 
not receiving over 2.5 times more fits for these CCU 
trends than the two control CCU trends. However, 
Belgaum received less fits to the 'estimated CCU trends' 
than the 'Control CCU trends' (Table 3). These results 
suggest that large increases in CCU most probably 
occurred during Avahan in all districts other than 
Belgaum and possibly Bangalore. 

Predicted female sex worker and client HIV prevalence 
trends over time 

Model projections based on the 'estimated CCU trends' 
suggest that FSW HIV prevalence has been decreasing in 
all districts since 2004-2005. Indeed, the FSW HIV 
prevalence may decline to low levels by 2015 (<2— 10% 
depending on district) without the need for any 
additional interventions (Fig. 1). However, prevalence 
peaked earlier in districts where Avahan was not the first 
intervention (Belgaum, Bellary) because condom use 
increased earlier than in other districts, and the Belgaum 
HIV epidemic was more mature at the start of Avahan 
than elsewhere (Supplementary figure S5, http:// 
links.lww.com/QAD/A319). Without condom use 
increasing after 2004 (Control 1 CCU trend), the 
FSW and client's HIV prevalence may have remained 
fairly stable or increased slightly for a few years in Mysore, 
Bellary and Bangalore; decreased in Belgaum; and 
increased in Shimoga (Fig. 1, Supplementary figure S6, 
http://links.lww.com/QAD/A319). 

HIV incidence and infections averted 

At the Rl IBBA, the modelled FSW HIV incidence was 
lowest in Shimoga (2.6 per 100 person-years) and highest 



inMysore (11.7 per 100 person-years), and reduced to 1.1 
and 3.3 per 100 person-years at R2/R3, respectively 
(Supplementary table S4, http://links.lww.com/QAD/ 
A3 19). Table 4 shows the predicted prevented fraction 
and number of HIV infections averted among FSWs and 
clients for each district over different time periods, 
estimated by assuming either counterfactual (control 
groups 1 or 2). Using Control 2, the median prevented 
fraction among FSWs (among clients) increased from 
16-20% (17-23%) in year 1 to more than 43% (>49%) 
over 5 years across all districts except Belgaum (Table 4B). 
The increase in condom use following Avahan is 
estimated to have prevented 36.7% (95% CrI 27.0— 
43.7%, Belgaum) to 55.4% (95% CrI 39.1-68.1%, 
Shimoga) of new HIV infections among FSWs from 
the start of Avahan to 201 1 , and between 46.0% (95% CrI 
34.2-54.0, Belgaum) and 59.0% (95% CrI 47.3-72.8, 
Bangalore) among clients (Table 4B). Due to differences 
in population sizes across districts, this translates to 
between 220 (95% CrI 120-406, Shimoga) and 1255 
(95% CrI 660-2530, Bangalore) new HIV infections 
being averted among FSWs by 2011. Two-fold to nine- 
fold more infections are averted in clients than FSWs 
(Table 4) in each district due to their larger population 
size. Lastly, syphilis treatment alone is estimated to prevent 
on average a maximum of 6% of new HIV infections over 
10 years (Table 4C). 

Discussion 



This is the first evaluation of the Avahan intervention 
among high-risk populations in Karnataka state using a 
preplanned evaluation framework [8] . Our results suggest 
Avahan has been effective, particularly in Bellary, 
Shimoga and Mysore. FSW HIV prevalence declined 
in all districts, with statistical significance in three districts. 
The decline in HIV prevalence may have been facilitated 
by the fast turnover of sex workers. STI trends were less 
consistent than HIV trends potentially reflecting variation 
in programme effectiveness, degree of implementation 
(e.g. syphilis screening took off in 2008, except for 
Mysore where it started earlier), and different and faster 
transmission dynamics or a lack of statistical power. 
However, delays in conducting the first IBBA mean that 
the prevalence of curable STIs (especially those with short 
duration such as chlamydia or gonorrhoea) may have 
already decreased appreciably before the Rl IBBA 
[21,22,35]. Thus, the observed reduction in STI may 
underestimate the true decrease following programme 
implementation. Interestingly, Mysore, which had the 
earliest Rl IBBA, noted the largest declines in all STIs 
between Rl and R2. Lastly, a recent analysis combining 
the five Karnataka districts showed significant declines in 
HIV and all STIs between Rl and R2 IBBA [7]. 

Our Bayesian modelling analysis strongly suggests that the 
observed declines in FSW HIV prevalence in Shimoga, 



Copyright © Lippincott Williams & Wilkins. Unauthorized reproduction of this article is prohibited. 



Evaluation of HIV prevention among high-risk groups Boily ef al. 1455 



> 
< 



o 

■D 



D 

U 

u 



E &) 



! o 9 



- < 2 

o CO 

cq "a 
~a — ^ 

CD < 

asm 



"D 5 "D 
CD CQ 

CO < 
= — CO 



ro OL - 



e a 



■o 

o 
o 



5 



cn O O 



, ro ? c c 0 



CN IT) 



cm <N 



^ ^ r_ 



^ o N ffl ° 

g rn CO - 



Ln * 

o- 



co \ 



<o o °"> 

U-l r- 



* CO 



ooooooooooooooo 
ooooooooooooooo 
ooooooooooooooo 
ooooooooooooooo 
ooooooooooooooo 

(NfNfNfNtNtNOOOi-'-r-lflinin 



< 


< 


< 


< 


< 


T3 




T3 




■c 


C 


C 


C 


c 


c 


CLJ 


QJ 


aj 


QJ 


aj 












D 


D 


Z3 


D 


D 


U 


U 


U 


U 


U 


U 


U 


u 


u 


U 


t— cm "a t— 


(N "D i— 


(n "a t— 


r-j "a t— 





o o 



QJ — — QJ„„ QJ. 



QJ . 



O O ' 



QJ 



_ .E" c c .e" c c M c c .e" c c .i 

oo-5;oot?500tsoo-KOOts 



QJ 



& o 
£ IE 



QJ t-O 

E .E 



e ; 

-s Si < 

&-.§>! ij. 



00 to OJ 



aj 



o cn 
En 



^ r XI 
CL t c 

E _oj qj 
o Q-i= 

U Q.— j 

% « U 



S £ 6 

QJ O <J 

E QJ DJD 

2 E .E 

re ~ 's£ 

Cl qj 3 



QJ 5 TO r 



QJ " 



o 

|| 0) TO • 

D c 

TO (J OJ 
> (J — 

°--0=> 

TO S U i 
* TO U 

H E £ 
re ™ c 

O < CN 

5 v>""e ' 

■Sec 

CD O - 

| fcU 

^ TO 

2^ Si 
^1 



TO r- 



E = 



cd j 



— . CO QJ 



~i re -= 

E - 
QJ .E c 

S « «3 

U J3J ^ 

li §_ 

+j re 

I'll 

S 2 > E 
re L. rt rs 

ill 

toU^ 

o a 2 

TO 

TO E D 

S »u 

H 1 

D£ -m .E 



— , « QJ 

qj ^ re 



DO qj qj 

C D0_o 

re ns c 

23 = 



re qj c 
E qj £ 

"Id 

DO TO -J 
"D < 7 

c cq a 

ro cq S 

I 

re 

03"°. oj 



< o ■ 



^ * QJ 

o 

HI3 

S o u 
-3! >- E 

TO TO 

E »m 

CD ^3 • • 
C 

^ CD — 
CD CJ U 



> c e a 

; n S _ 

I "D CD 

j ' ' CD 

03 S E 

£9 to „ 
ro 



t > m cd 



ro ro 

— -a cr 

CD CD CD 

-C CD "O 

^ c ro 

7= CD 

c 15 0 

S £ s 

ro °- ez 

~ « CD 

S ' ro « 

Si > H 
£ to 

S 3! o 



o c 

CD CD 
. CD - 

3 s 
u a " 

U Jj : 



'cd to 

S ■ - ^ 

S"S CD 

Si Si o 
ro 00"^ 
cd ro go 

CD "D C 
C CD TO 

•-O.CQ 

i8« 



Copyright © Lippincott Williams & Wilkins. Unauthorized reproduction of this article is prohibited. 



1456 AIDS 2013, Vol 27 No 9 



(a) Mysore 



■ Estimated CCU trend 

■ Control 1 CCU trend 

■ Control 2 CCU trend 



0 
CO 
> 
CD 

Q. 

> 
X 

05 




1985 1990 

(b) Belgaum 



1995 



2000 



2005 



2010 2015 



• Estimated CCU trend 

■ Control 1 CCU trend 

■ Control 2 CCU trend 



I 

> 
X 



05 

LL 




1985 



1990 



1995 



2000 



2005 



2010 2015 



(c) Bellary 









O 


CD 




O 




c 




CD 






o 




CO 


<D 




Q. 




> 


o 


X 


C\J 


<: 




CO 




LL 


o 



. Estimated CCU trend 
. Control 1 CCU trend 
■ Control 2 CCU trend 




1985 1990 1995 

(d) Bangalore urban 



2000 



2005 



2010 



2015 



_CD 
CCS 
> 
CD 

> 
X 



. Estimated CCU trend 
. Control 1 CCU trend 
. Control 2 CCU trend 




1985 



1990 



1995 



2000 



2005 



2010 



2015 




Fig. 1. Female sex worker HIV prevalence over time simulated under the three different condom hypothesis trends. Predicted 
female sex worker HIV prevalence over time simulated under the 'estimated consistent condom use (CCU) trends' in presence of 
the intervention and the two alternative CCU trends (Control 1 and Control 2 CCU trends) to simulate HIV prevalence in absence of 
the intervention in (a) Mysore, (b) Belgaum, (c) Bellary, (d) Bangalore Urban and (e) Shimoga districts. Shown on the graphs are the 
mean (dark lines black, blue and red) and the 75% credibility interval (Crl; shaded area) for each hypothesis. The paler grey area 
represents the 95% Crl. Also shown is the available Integrated Behavioural and Biological Assessment (IBBA) survey prevalence 
data (mean and 95% confidence intervals, CIs). As the model was fitted to the adjusted HIV prevalence in R2 and/or R3, the HIV 
prediction may not fall exactly within the 95% CI of the HIV prevalence data. The transient increase in prevalence with Control 1 
and 2 CCU trends in (a) and (c) are due to the increase in number of clients of FSWs after R1 . 
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Table 4. Predicted number (No) and cumulative fraction (prevented fraction) 3 (95% credibility interval) of new HIV infections averted due to 
the increase in condom use and the improved syphilis treatment since the start of Avahan comparing model predictions with the 'estimated 
consistent condom use trends' and Control group 1 or 2. d,e 



A. Using control group 1 b 




Over 1 st year 


Over 3 years 


Over 5 years 


To end of 201 1 


To end of 2014 


Mysore — FSW 


PF (%) 


17.1 (11.7-24.3) 


35.9 (30.2-41 .3) 


47.4 (41.5-54.3) 


54.1 (47.7-61.2) 


61.1 (54.4-68.2) 


Shimoga —FSW 


PF (%) 


18.6 (12.3-27.4) 


40.7 (29.6-53.2) 


59.1 (47.5-69.0) 


66.9 (56.0-75.6) 


76.8 (67.7-84.0) 


Belgaum —FSW 


PF (%) 


19.9 (15.1-24.1) 


36.4 (28.6-42.1) 


43.9 (36.4-49.1) 


47.8 (40.9-52.8) 


53.5 (46.7-58.0) 


Bellary -FSW 


PF (%) 


22.9 (14.4-29.2) 


46.1 (31.9-56.5) 


60.9 (48.0-71 .1 ) 


65.4 (52.4-75.3) 


71 .4 (59.3-80.5) 


Bangalore Urban —FSW 


PF (%) 


21.2 (15.0-27.9) 


42.2 (32.2-52.7) 


54.2 (43.6-64.4) 


57.4 (46.4-67.2) 


67.3 (56.5-76.4) 


B. Using control group 2 b 




Over 1 st year 


Over 3 years 


Over 5 years 


To end of 201 1 


To end of Avahan (end of 2014) 


Mysore -FSW 


PF (%) 


15.6 (10.4-22.5) 


31.7 (26.4-39.1) 


42.6 (36.5-50.8) 


49.2 (42.4-57.6) 


56.2 (48.7-65.0) 




No 


30 (19-47) 


187 (126-272) 


411 (273-592) 


633 (409-926) 


1004 (616-1476) 


Clients 


PF (%) 


17.0 (11.1-24.3) 


36.3 (30.4-46.2) 


49.0 (42.3-57.8) 


55.5 (48.1-64.3) 


62.2 (54.1-70.9) 




No 


70 (34-133) 


475 (236-769) 


1033 (502-1644) 


1512 (706-2511) 


2221 (994-3851) 


Shimoga -FSW 


PF (%) 


15.8 (9.2-25.5) 


33.7 (21.7-46.4) 


48.3 (33.6-61.1) 


55.4 (39.1-68.1) 


65.7 (48.0-77.8) 




No 


9 (5-18) 


60 (32-99) 


148 (80-257) 


220 (120-406) 


400 (218-769) 


Clients 


PF (%) 


17.0 (10.4-26.6) 


35.8 (23.3-47.9) 


51.6 (36.9-63.5) 


58.9 (43.6-71.0) 


68.9 (53.2-80.6) 




No 


58 (25-110) 


356 (158-671) 


835 (405-1566) 


1228 (594-2506) 


2189 (1004-4619) 


Belgaum -FSW 


PF (%) 


17.8 (12.4-21.3) 


29.8 (22.1-36.2) 


34.1 (25.1-40.6)% 


36.7 (27.0-43.7) 


40.5 (29.4-48.3) 




No 


23 (17-32) 


1 15 (86-1 56) 


206 (142-273) 


276 (183-365) 


393 (248-524) 


Clients 


PF (%) 


21.5 (15.9-27.2) 


38.9 (29.3-47.9) 


43.4 (32.4-52.1) 


46.0 (34.2-54.0) 


49.4 (36.1-57.1) 




No 


209 (138-352) 


1047 (710-1697) 


1 791 (1 129-2863) 


2326 (1444-3678) 


3216 (1989-4986) 


Bellary -FSW 


PF (%) 


19.5 (11.8-26.7) 


38.1 (24.1-49.1) 


49.0 (32.8-62.1) 


52.7 (35.6-66.4) 


57.9 (39.3-72.2) 




No 


27 (13-51) 


161 (75-285) 


343 (167-615) 


444 (214-800) 


633 (288-1098) 


Clients 


PF (%) 


21 .7 (12.9-28.4) 


42.7 (27.7-52.3) 


54.5 (38.5-66.4) 


58.6 (41 .4-70.9) 


63.8 (44.6-76.1) 




No 


257 (115-532) 


1490 (605-31 70) 


3152 (1387-7465) 


4106 (1748-9405) 


5497 (2255-11 710) 


Bangalore Urban -FSW 


PF (%) 


19.5 (13.6-27.8) 


37.5 (27.6-52.5) 


47.9 (37.5-63.5) 


50.6 (40.1-66.2) 


60.1 (47.2-75.1) 




No 


93 (55-158) 


536 (311-948) 


1055 (574-2092) 


1255 (660-2530) 


2220 (1024-4822) 


Clients 


PF (%) 


23.2 (16.3-31.7) 


44.6 (32.9-59.1) 


56.1 (44.5-70.3) 


59.0 (47.3-72.8) 


68.0 (56.8-81.0) 




No 


295 (139-558) 


1744 (722-3465) 


3476 (1325-7506) 


4104 (1517-8985) 


6750 (2314-16660) 


C. Tp treatment only c 




Over 1 st year 


Over 3 years 


Over 5 years 


To end of 201 1 


To end of Avahan (end of 2014) 


Mysore -FSW 


PF (%) 


3.1 (1.0-7.5) 


2.3 (-0.2-7.5) 


2.4 (-0.2-8.2) 


2.6 (-0.3-8.6) 


2.7 (-0.5-9.4) 


Shimoga -FSW 


PF (%) 


3.7 (0.3-7.8) 


4.8 (0.3-10.1) 


5.2 (0.3-11.4)% 


5.3 (0.3-12.3)% 


5.6 (0.3-13.6) 


Belgaum -FSW 


PF (%) 


3.0 (1.0-6.3) 


2.8 (0.4-6.4) 


2.7 (0.1-7.2) 


2.8 (0.0-7.8) 


3.1 (-0.1-8.8) 


Bellary -FSW 


PF (%) 


5.0 (2.1-9.1) 


5.6 (1.8-12.0) 


5.3 (1.4-13.2) 


5.4 (1.3-13.4) 


5.4 (1.4-13.2) 


Bangalore Urban -FSW 


PF (%) 


3.5 (0.9-6.8) 


3.9 (0.4-7.7) 


4.0 (0.3-8.6) 


4.0 (0.2-8.9) 


4.3 (0.1-10.0) 



CCU, consistent condom use; FSW, female sex worker; PF, prevented fraction. 

Prevented fraction, [number of new infections occurring between the start of Avahan (T1 ) and time T2 in absence of the intervention (h) minus the number of new 
infections between T1-T2 in presence of the intervention (l 0 )] divided by h. 
b Reflects the Impact of Syphilis (Tp) treatment and condom use. 

c Reflectsthe Impact of Syphilis (Tp) treatment only - assumes same condom use as the 'estimated CCU trends' but no increase in Tp treatment due to Avahan - note that 
negative PF values are possible because in some instances, the more rapid syphilis treatment makes the latent people susceptible again that can increase syphilis 
incidence and prevalence over time and increase HIV incidence through the cofactor effect for HIV. 

d AII these projections assume constant condom use after the last round of Integrated Behavioural and Biological Assessment (IBBA) surveys (R2 in Shimoga and 
Bangalore, R3 in Mysore, Belgaum and Bellary). 

Impact estimates and 95% credibility intervals (95% Crl) were generated by using the projections from each posterior parameter set weighted by the likelihood of the 
model projection as explained in methods and supplementary material, Methods C, http://links.lww.com/QAD/A319. 



Mysore and Bellary are consistent with large self-reported 
increases in condom use during Avahan ('estimated CCU 
trends'). The evidence was moderate for Bangalore and 
weaker for Belgaum. Thus, CCU increases following 
Avahan are likely responsible for most of the decline in 
FSW HIV prevalence in all districts, except Belgaum. In 
Belgaum, HIV prevalence may have declined without 
condom use increasing during Avahan, partly due to the 
natural dynamics of the epidemic and pre-Avahan 
increases in condom use because Avahan was not the 
first intervention. This highlights how models can help 
interpret HIV time trends more cautiously when 
evaluating an intervention's impact. 

Our modelling analysis highlights the time required to 
substantially avert infections among HRGs. Across the 
four districts with stronger evidence (excludes Belgaum) 
for increasing condom use during Avahan and using a 
conservative counterfactual (Control 2), we estimated 



that 36—68% of new HIV infections were prevented 
among FSWs over 7 years (until 2011) compared to 
9-28% and 22—53 over 1 and 3 years, respectively. 
Interestingly, up to 20% larger fraction of infections were 
averted among clients. 

Given the 'real life' evaluation challenges faced (no 
randomized control groups, delayed baseline measure- 
ments), our results have limitations, which have been 
partly circumvented through our analytic approach. First, 
our 'estimated CCU trends' in presence of intervention 
were based on reconstructed CCU trends from self- 
reported behavioural data, which were liable to recall and 
social desirability biases [9]. However, our Bayesian 
modelling analysis showed that these 'estimated CCU 
trends' are generally more consistent with HIV preva- 
lence trends than assuming no or slow increases in CCU, 
suggesting that these trends have some basis. Despite the 
modelling analysis suggesting that these increases in 
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condom use partially caused the observed decline in HIV 
prevalence, it cannot directly attribute the CCU 
changes to Avahan without additional information 
(Supplementary material, Background section, http:// 
links.lww.com/QAD/A319). Complementary process 
monitoring data collected as part of Avahan indicate 
that the intervention has been rapidly implemented and 
scaled up in Karnataka (Supplementary material, Box li, 
http://links.lww.com/QAD/A319) [10,11]. Consistent 
with our 'estimated CCU trends', additional evidence 
suggests that Avahan may have accounted for 88% of the 
abrupt increase in condom availability following the start 
of Avahan in these Karnataka districts (Supplementary 
material, Box lii, http://links.lww.com/QAD/A319) 
[36]. Importantly, analyses of Rl and R2 FSW 
behavioural data from Karnataka found a dose— response 
relationship between intervention exposure and self- 
reported condom use in commercial partnerships 
(Supplementary figure S4, http://links.lww.com/ 
QAD/A319) [37J. In this dose— response analysis, the 
levels of CCU with occasional clients reported by those 
not in contact with the Avahan intervention were 
compatible with the CCU levels obtained from the 
average 'estimated CCU trends' around the time that 
Avahan started (Supplementary figure S3-S4, http:// 
links.lww.com/QAD/A319). Together, these results 
suggest that much of the CCU increases can be attributed 
to Avahan, especially in districts where Avahan has 
remained the only HRG programme (Mysore, Shimoga). 

To minimize model misspecification, the model structure 
was carefully developed to reflect important sources of 
observed behavioural heterogeneity, based on analyses of 
the key determinants of HIV infection in the IBBA 
surveys. Nevertheless, as with any modelling analysis, 
simplifying assumptions were made. As no data on sexual 
mixing were available for commercial sex work partner- 
ships, we assumed random mixing between clients and 
FSWs. Although sexual mixing influences the spread of 
infection, its effect on intervention impact projections is 
more modest when models of high-risk populations are 
fitted to detailed HIV prevalence data [38] . In the absence 
of HIV incidence data, we fitted the model to FSW HIV 
prevalence data, and estimated infections averted using 
modelling. The modelled HIV incidence at the time of 
the first IBBA ranged between 2 and 12 per 100 person- 
years across districts, reflecting different intervention 
histories and condom use levels prior to Avahan, which is 
not directly comparable but is consistent with the only 
available HIV incidence estimate from Pune between 
1993 and 2007 [39]. Our Bayesian framework allows to 
reflect parameter uncertainty in our impact estimates — 
using only one parameter set could have produced overly 
optimistic or pessimistic estimates. Finally, in validation 
experiments, our hypothesis testing methods showed 
good discriminating power to identify or reject increasing 
CCU trends when it is true or false, respectively, 
especially for districts with four IBBA rounds or where 



the ratio of fits exceeds one appreciably (Supplementary 
table S5, http://links.lww.com/QAD/A319). 

Unfortunately, comparable data on HIV/STI prevalence 
and behaviour among FSWs over time were only available 
from intervention districts [3,8,14,15]. Given the lack of 
control group to indicate how condom use would have 
evolved without Avahan, we based our impact estimates 
on a conservative counterfactual, which assumed condom 
use would have increased at preintervention rates without 
Avahan. This is more conservative than assuming a stable 
condom use counterfactual as typically done. Despite the 
use of mapping and cluster random sampling, changes in 
the composition of the FSW population may have 
influenced the decline in HIV prevalence. However, we 
accounted for potential differences in two ways: by 
allowing the parameters for the frequency of commercial 
sex acts to change over time in the modelling analysis, and 
by accounting for other changes in key confounding 
factors by using adjusted HIV prevalence at the fitting 
stage. There was no evidence for a large reduction in FSW 
duration across IBBA rounds. We did not model ART 
because coverage was low during Avahan and unlikely to 
have influenced results. 

Our analysis has several implications for future evaluations 
of HIV intervention programmes. Here, the Bayesian 
modelling approach was crucial for the success of the 
evaluation. Otherwise, the evaluation would have relied 
on trends in risk behaviour and HIV/STI prevalence from 
delayed IBBA surveys only. Our modelling analysis 
helped interpret trends and estimate impact while 
accounting for the natural dynamics of infection, past 
history of interventions, district-specific epidemiological 
contexts and uncertainty in the 'estimated CCU trends'. 
The use of models for hypothesis testing provided a 
mechanism to grade the strength of evidence for 
intervention impact. 

Our results have important policy implications. The 
empirical and modelling results provide plausible evidence 
that condom use during commercial sex has increased 
during Avahan, and contributed to reducing HIV 
transmission between FSWs and clients in Karnataka. 
These assertions are strengthened by strong evidence for 
the successful implementation of the Avahan programme, 
and that increases in condom use among FSWs are 
associated with programme exposure. Our study supports 
the notion that HIV prevention programmes targeted at 
HRGs are feasible and can have considerable impact. 
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